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Abstract 
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Chern-Simons-Higgs vortex equations over a doubly periodic domain when the gauge 
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1 Introduction 

Let K = ( Kij ) be the Cartan matrix of a semi-simple Lie algebra L. Recall that the Toda 
system is a system of nonlinear elliptic equations over M 2 , of exponential nonlinearities, of 
the form 

r 

A Uj = i=l,...,r, (1.1) 

3 = 1 

where r is the rank of L. This system is known to be integrable in general [^111221131115411421147] 
and arises in the study of non-Abelian monopoles [22j,;37],[53j and nonrelativistic Chern- 
Simons-Higgs vortices [17314161 [25] . Interestingly, when r — 1, that is, when the Cartan 
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subalgebra of L is Abelian such as when L is the Lie algebra of SU( 2), (II. lj) reduces to the 
classical Liouvillc equation [38] 

A u = —Xe u , (1.2) 

whose solutions may be constructed by all the integration methods known, such as separation 
of variables, inverse scattering, the Backlund transformation, etc., and is often used as an 
illustrative example. It is well known that the Liouville equation and its extensions arise 
also in differential geometry [2ll5[l29[l3Ci] and have been the focus of various studies on their 
analytic aspects |3ll6L fT2lf35] . More recently, some extensive work on (11.11) has been carried 
out as well aimed at the classification of solutions [28jl36j, understanding its fine analytic 
structures [39l ITT] . and establishment of bubbling behavior of solutions [2711371145 l| . 

On the other hand, however, when one considers relativistic Chern-Simons-Higgs vortices 
[2,26, 56], the governing system of equations [I^l - [l6l l56] is 

r r r 

Au * = A^^AV^“ J e u *i — 1,... ,r, (1.3) 

j =1 k =1 j =1 

which deviates from (11.11) significantly in that, even in the scalar case where r — 1, the 
system is nonintegrable [48] . Thus, it is perceivable that the analytic structure of (11.31) 
would be much more complicated than that of (11.11) . Due to its applicability in anyon 
physics PM233E3] and challenging mathematical content, it will be desirable to develop 
an existence theory for the relativistic vortex equations (11.31) , in which the sources terms 
resulting from the presence of vortices are temporarily neglected in order to facilitate our 
discussion. 

In [55], a systematic study is conducted to establish an existence theorem for the so-called 
topological solutions, realizing the spontaneous symmetry breaking or the celebrated Higgs 
mechanism of the model, over the full plane K 2 , to the equations (11.31) . where the Cartan ma¬ 
trix is of a general form. Another type of solutions of great interest are called the Abrikosov 
vortices p] or vortex condensates for which the equations are to be solved over a doubly- 
periodic lattice domain. In field-theoretic formalism, such a structure is realized by imposing 
the ’t Hooft [50| boundary condition on gauge and matter fields [I9]l52j|56]. Surprisingly, 
although the system is now considered over a compact domain, double periodicity greatly 
complicates the problem as evidenced [46] already in the nonrclativistic Liouvillc equation 
situation where one needs to use the Weierstrass elliptic functions as building blocks for 
solutions. In the relativistic situation, although there is a variational principle, the action 
functional is not bounded from below and one has to consider a constrained minimization 
problem. However, the presence of the constraints leads to a Lagrange multiplier issue so 
that it prevents one from recovering the original equations of motion. In order to overcome 
this difficulty, one has to consider an inequality-constrained problem instead and bypass the 
Lagrange multiplier problem with achieving an interior minimum. As a consequence, the 
progress in developing an existence theory for the doubly periodic solutions of the relativis¬ 
tic Chern-Simons-Higgs vortex equations has been slow and sporadic. Specifically, in [3], 
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an inequality-constrained minimization method was first used to establish the existence of 
solutions for the scalar case of (11.31) . namely r = 1 or G = £77(2), and, in pH], the method 
in [3] was remarkably extended and refined to tackle the first non-scalar case of (1 1.31) . namely 
r = 2 or G = SU( 3), in which there are two inequality constraints characterizing the solv¬ 
ability of two quadratic constraints. The study on the next important non-scalar case of 
o. that is, r = 3 or G = 577(4), had been unsuccessful due to the difficulty in resolving 
more than two quadratic constraints simultaneously until very recently a new idea based on 
an implicit-function theorem argument was implemented to resolve three coupled quadratic 
constraints, which enables the establishment of an existence theorem for G = £77(4) in [231. 

The purpose of the present paper is to establish an existence theorem for the doubly 
periodic solutions of (II. 3ft for the most general situation, G = SU(N) (N > 2) or r = 
N — 1 > 1, adapting the implicit function method initiated in [23J in handling multiple 
quadratic constraints. Instead of a ‘squeeze-to-the-middle’ approach in [23] for resolving the 
constraints, however, we use here an ordered iterative scheme which is effective and easier to 
implement in the general situation. In the next section, we state our main existence theorem. 
In the section after, we prove the theorem. In the last section we end the paper with some 
remarks. 

Note added upon acceptance. This work grew out of an earlier version of [23] in which 
only the £77(4) problem was resolved. Since the submission of the present work, further 
development has been achieved to resolve the general situation when G is a simple Lie group 
by using a degree-theory method formalism, which greatly expands [23] into its current 
updated version seen. 


2 Existence theorem 


With the source terms in the presence of multiply distributed vortices, the non-Abclian 
relativistic Chern-Simons-Higgs equations are HMM 


( r r r \ Ni 

Y Y Kk j K ji e Uj e Uk - Y K ^ UJ )+^Y 6 ph i = h---,r, ( 2 . 1 ) 

3 = 1 fc=l 3 = 1 / 3 = 1 

where S p denotes the Dirac measure concentrated at the point p, A > 0 is a coupling constant, 
and the equations are considered over a doubly periodic domain resembling a lattice cell 
housing a distribution of dually charged vortices located at p l3 , j = 1 ,..., = 1 ,,r. 

We focus on the system (12.ip with G = SU{n + l),n > 2. The Cartan matrix K for 
SU(n + 1) is an n x n matrix given by 


/ 2 -1 0 ••• 0 \ 

-1 2 -1 ... 0 


K = 


0 

\0 


-1 2 -1 
0 - 12 / 


( 2 . 2 ) 
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It is easy to check that K 1 is symmetric with entries given by 

(K-% = i<j,i = l,...,n. (2.3) 

n + 1 

Then after the translation 

Ui y Ui T In/?,, with A,; = =---, i = (2.4) 

3 =1 

the system (12.11) can be rewritten as 


Ni 


+47ry^.(x), z = l,...,n, (2.5) 

v j=i fe=i i=i / t=i 


or in a vector form, 


where 


Au = AATJ/T(U - 1) + 4tts, 


( 2 . 6 ) 


/ n —{n — 1) 
-f 2{n — 1) 


K = I\R = 


0 

0 


0 

3( n —2) 


■ ( ^ 1)( " +2 ~^ i(ra + l-i) -. (i+1 K"-0 . 


0 \ 
0 


2(n — 1) 


, (2-7) 


-(n — 1) n y 


A = diag{i?i,. . . ,i? n } , 1 = (1,..., l) r , 


u= (u 1 ,...,u n ) T , U = diag{e Ml ,... ,e Un } , U = (e Ul ,..., e^ 


JVi 


N n 


S _ I > • • • ! V. 


( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 


, S=1 


s=1 


By the dehnition of K given in (12.71) . we obtain the following simple facts 
R- 1 = R-'K- 1 , K~ l l = 1, AT" 1 ! = Rl, 


( 2 . 11 ) 


which will be repeatedly used later in this paper. 

We are interested in the existence of solutions of (12.51) or (12.61) over a doubly periodic 
domain hi. Our main result reads as follows. 
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Theorem 2.1 Consider the nonlinear elliptic system (j2.5[) or (I2.6[) over a doubly periodic 
domain 12 in M 2 . For any given points pn,... ,PiN t £ 12 (z = 1,... n), which need not to be 
distinct, the following conclusions hold. 

(i) (Necessary condition for existence) If 


A < A 0 = 


EE ( k -% n , 

l07T i=1 j =1 


VI 


EEin 

i=lj=l 


( 2 . 12 ) 


there is no solution to the system, where the entries of K l are given by (12.3D . In other 
words, a solution can exist only when A is larger than the right-hand side of (12.121) . 

(ii) (Sufficient condition for existence) There exists some Ai > Ao such that when A > Ai 
the system admits a solution over 12. 

(iii) (Asymptotic behavior) The solution (u \,... ,u n ) obtained above satisfies 


lim / (e Ui — l) 2 dx = 0, i = 1,..., n. (2.13) 

A->oo J n 

(iv) (Quantized integrals) If {u\,... ,u n ) is a solution then there hold the quantized inte¬ 
grals 



(2.14) 


Note that when the vortex numbers Ni(i = 1,,n) are the same, say m, our necessary 
condition (12.12[) reduces to that for the U( 1) case as in [4] as follows 


A < 


IGnrri 

W 


(2.15) 


which is a comfort and also surprising since we are now considering a non-Abelian and 
non-scalar situation. 


3 Proof of theorem 

In this section we apply a constrained minimization procedure developed in [0], which was 
later modified in [44] . to establish the existence of doubly periodic solutions to (12.5p (12.6[) 
when n — 2. We carry out our proof in several steps. First we show that the condition 
(12.12)1 implies nonexistence of solutions as stated and we explore a variational structure of 
our equations. Next we show how to resolve multiple constraints using an iterative scheme 
and an implicit function argument. We then conduct a constrained minimization procedure 
and show that there is a solution to the minimization problem. In the subsequent subsection, 
we establish some suitable estimates which ensures that the minimum point obtained must 
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be an interior minimum, thus ruling out the Lagrange multiplier issue. In the last two 
subsections, we show that the interior minimum point obtained is a classical solution of the 
original vortex equations and we then establish the stated asymptotic behavior of solutions 
as A —> oo and quantized integrals. 


3.1 Necessary condition and variational structure 

Let u® be the solution of the following problem (see 0) 


Ni 


Au°i = 4n S P is 


tqdx = 0, i — 1,..., n, 


S=1 


M ’ ^ 

and Ui — + Vi, i — 1,... ,n. Let us introduce the notation (n- vectors) 

v = (ui,... ,v n ) T , N = (iVi,...,A^ n ) T , 0 = (0,..., 0) T . 

Then we reformulate the system (12.5p or (I2.6[) as 


\j =1 fe=l 


lJ e u °i +Vj e u ° +Vk -J2^ij eU ° j+Vj I + 


3 = 1 


VI 


* = l, 


(3.1) 


(3.2) 


(3.3) 


(3.4) 


(3.5) 


We hrst present a necessary condition for the existence of solution to (13. 3 j) or (13.41) . For 
any solution v of (13.4[) . taking integration over hi, we obtain the natural constraint 


or 


47T 

Av = \KXJK(XJ - 1) + —N, 

with the understanding that 

U = diag { e u ° +vi ,..., e u ° +Vn } , U = (e^ 1 ,..., e <+ ^) T . 


f ~ ~ in 

/ K\JK(XJ - l)dx + —N = 0. 
In * 


(3.6) 


Multiplying both sides of (13. 6 p by K l , we have 

C Ajr 

/ RXJKR(XJ - l)dx + —iW 1 N = 0. (3.7) 

Jn A 

Then noting (13.7p . the positive definiteness of K, and the fact that K~ l 1 = Rl, we have 


r 4-7T 

0 = (RXJ) T KR{XJ-l)dx + —l T K- 1 N 
Jn A 


( 

\ 1 1 

V ( 

\ i 1 

[n 


) K \ R 

[ U _ 2} 


da; - ^(R1) T K(R1) + y VK^N 


> -^1 T K- 1 1 + —1 T K~ 1 N, 
4 A 


(3.8) 
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which implies 


(3.9) 


A > 


167rl I\ x N 

|fi|l T /l _1 1 


That is, (13. 9 p spells out a necessary condition for the existence of solutions to (13.4p . Hence 
the first conclusion of Theorem 12.11 follows. 

Now we aim to find a variational principle for the equations (13.4p . After a simple com¬ 
putation, we see that the matrix K admits a decomposition 


K = PS, 


(3.10) 


with 


A 

= diag |Pi, 

p \ 

' • ’ i ± n C ) 

p — " 


l ~ i{n + l- 

-i)' 


f 

n 

T— 1 

1 

T 

0 



~( n ~ 1) 

4(n-l) 2 

cT 

1 

7 



n 

n 


s = 



— Sii-l 

~Su 






3(n—l)(n—2) 

4(n—l) 2 





n 

n 


V 

0 


0 

i 

's' 

1 

H 

tridiagonal matrix with 





' Su-1 
Jii 

k Sii+1 

_ (i— l)i(n+2— i) 




i 2 {n+ 1- 

2 n ’ 

-h 2 

= 3, 


1 

n ’ 

_ i)(n—i) 

~ 2 n > 


0 \ 


n 


Then we may rewrite (13,4h equivalently as 


A Mv = AUA(U — 1) + 


M’ 


where 


(3.11) 


(3.12) 


(3.13) 


(3.14) 


and 


M = P- X K~ X = -K~\ P~ l = diag <1 P{ 


-1 D-l 


n 


3-1 


P~ l , 

? * * * > J n [ 5 


A- 1 = 


i[n + 1 — i) 


n 


, i = l,...,n, 


b = {b 1 ,...,b n ) T = 4nMN. 


(3.15) 

(3.16) 

(3.17) 
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By the definition (13.171) for b, we easily find that bi > 0, i — 1,..., n. 

We use W 1,2 (r2) to denote the Sobolev space of scalar-valued or vector-valued fl-periodic 
L 2 functions with their derivatives also belonging to L 2 (Q). 

ft may be examined that the equations (I3.14j) are the Euler-Larange equations of the 
action functional 

J(v) = ]- V' f diV T Mdivdx + ^ f (U — 1) T £(U — l)dx + —f b T vdx, (3.18) 

2 . =1 ,/n 2 J n J n 

where we use the notation (I3.12p . (13.15[) and (13.17[) throughout this paper. 

In the following subsections we will use a constrained minimization approach to find the 
critical points of the functional I. 


3.2 Multiple constraints 


To start our constrained minimization process, we need to find some suitable constraints 
subject to which the functional / will be minimized. 

Note the space iy 1,2 (f2) can be decomposed as 


W 1,2 (Q) 


where 


W 1>2 (Q) 


w 


e W h2 (Q) 




is a closed subspace of VE 1,2 (f2). 

Then, for Vi 6 lE 1,2 (f2), we have the decomposition 


Vi = Ci + Wi, q e t, Wi G hE 1,2 (t2), i — 


(3.19) 


(3.20) 


(3.21) 


If v G W 1,2 (t2) satisfies the constraint (13.6p . which is equivalent to 


U5(U- l)dx+ - = 0, 

A 


then by the decomposition (13.211) with w = (nq,..., w n ) T , we obtain 


(3.22) 


e 2ci a n - e cl Pi(w; e C2 ) + —h = Q, 


bi 


n\ 


e 2ci au - e Ci Pj(w;e Ci " 1 ,e Ci+1 ) + — 


nbi 


i 2 (n + 1 — i) 2 A 


= 0, % — 2,..., n — 1, 


e‘ n a nn - e Cn P n (w; e 0 "- 1 ) + — = 0 

nX 


(3.23) 

(3.24) 

(3.25) 
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where and in the sequel the notation 


Pi(w;e C2 ) =-b 

n 


(n - l)a 12 


vC2 


n 


Pi(w; e Ci_1 , e Ci+1 ) = 


a, (i- l)(n + 2-t)g ii _i cC ._ i 

i(n + 1 — i) 2 i(n + 1 — i) 


(i + l)(n - i)au+i 
2i(n + 1 — i) 


e Ci+1 , i = 2,..., n — 1, 


P n (w;e c ’ 1 - 1 ) = 


_ a n . 1 Cn _j 

r e , 

n n 


and the definition 


a* = ai(wi) = / e Ui+Wi dx , 
Jn 


Oy = «h) = / e“ i+ ^ +Wi+Wj 'dx, i, j = 1,..., n, 


(3.26) 

(3.27) 

(3.28) 

(3.29) 

(3.30) 


are used. 

Then, we see that, for any w € M /1,2 (f2), the equations fl3.23ll - fl3.25ll are solvable in 


only if 


Pit 

w;e C2 ) 

> 

f(w;e c< - 1 

> e Ci+1 ) 

> 

^n( W 

; e Cn_1 ) 

> 

can be ensured 

by 


a 2 

> 


c — (ci,..., c n ) , 


46ian 


nA 


Anbida 


i 2 (n + 1 — i) 2 A 
1 b n a n n 


, i — 2,... ,n — 1, 


nX 


Anb l a l 

A 


or expressed explicitly as 


e .;«. d xl > hb 


i = 1,... ,n, 


e 2u ° +2wi dx, i — 1,... ,n. 


(3.31) 

(3.32) 

(3.33) 

(3.34) 
;s 

(3.35) 


Now we define the admissible set 


A = < w 


w 


G W 1,2 (yt) satisfies (13.3611 j>. 


(3.36) 


(3.37) 
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Thus, for any w e A, we can find a solution of the equations (]3.23j) - (j3.25j) with respect 
to (c 1} , c n ) by solving the following coupled equations 


e Cl = 


e Ci = 


Pi( w;e°) + ypf( w; e o)-^T 


2a 


n 


/i(e C2 ), 

P,(w; e Ci 1 , e Ci+1 ) + y^P^w; e^- 1 , eCi+1 ) — 


4nb^ a 


z 2 (n+l—i) 2 A 


2 ( 2 7 ;- 


e = 


/i(e c<_1 , e Ci+1 ), z = 2,..., n — 1, 
P n (w;e c "- 1 ) + ^/P n 2 (w;e c -i) - 


2(i, 


= fn(e c - 1 ). 


(3.38) 


(3.39) 


(3.40) 


For any w e *4, to solve the system of the equations (j3.38[) (I3.40p with respect to 
(ci,..., c n ), we will convert the system into a single equation. 

For n > 2, we need to find a solution of the system 


ti~fi(t 2 ) = 0, (3.41) 

U ~ = 0, i — 2,... ,n — 1, (3.42) 

tn ~ fnit n - 1) = 0. (3.43) 


Proposition 3.1 For any w e A and n > 2, the system (13. 41 p - (13.431) admits a unique 
solution in (0, oo) ?l . 


Proof. When n = 2, 3, the system can be transformed into a single equation directly, 
which was solved in 031 and [23j, respectively. However, for n > 4, it difficult to reduce 
the system (j3.4ip ~ ()3.43p into a single equation directly. Here we use the implicit function 
theorem to overcome this difficulty. 

For n > 4, using the first two equations of the system (I3.4ip - (l3.43p . we have the relation 
F 2 (t 2 ,t 3 ) = t 2 ~ f 2 (fi{t 2 ),t 3 ) = 0, t 2 ,t 3 > 0. (3.44) 

We first show that the relation (I3.44p may uniquely determine an implicit function 

t'2 = <72(P) >0, t 3 > 0. 

By the expression (13.38p (13.40p . we have 

dfi{t 2 ) = 1 ( (n-l)oi 2 | (n ~^ )ai2 ^i(w;t 2 ) \ 

dt 2 2 a n | n ^(w; t 2 ) - ^ f 

\MV; y - ^ 


(3.45) 


(3.46) 
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Similarly, we obtain 


df2(h,t 3 ) _ 4^=T)/2(^l,^3) 


dtl ^(wiMs)-,^' 

Then, using (13,46j) 113.47H and the constraints (13.3511 , we have 


(3.47) 


dF 2 (t 2 ,t 3 ) = 1 _df 2 (t 1 ,t 3 ) 


dt 2 


dti 


= 1 


> 1 - 


d/i(t 2 ) 

ti=fi(t2) d t 2 

la 2 12 f 1 (t 2 )f2(fi(t 2 ),t 3 ) 


\JpK W ;«2) - 

\al 2 fi{t 2 )f 2 (fi{t 2 ),t 3 ) 

(n— l)ai2 + nai2 




= ! _ MMiM = hfaM = o, t 2 ,t 3 > o. 

to to 


(3.48) 


Therefore, from (13.4811 and the implicit function theorem we see that there exists a unique 
implicit function 


h = g 2 (t 3 ) >0, t 3 > 0, such that F 2 {g 2 (t 3 ), t 3 ) =0, t 3 > 0. 


(3.49) 


Similarly, combining (j3.49j) and the third (i = 3) equation of the system (13.4111 - 03.4311 
gives the relation 


F 3 (t 3 ,t 4 ) = t 3 — f 3 (g 2 (t 3 ),t 4 ) — 0 , t 3l t 4 > 0 , 
which also determines a unique implicit function 

t3 = .9.3(( 4 ) >0, t 4 > 0, 

such that 


(3.50) 


(3.51) 


^ 3 ( 93 ( 14 ) 44 ) = 0, 1 4 > 0. (3.52) 

Repeating the above procedure, we obtain that there exists a family of uniquely deter¬ 
mined implicit functions 

ti = gi(t i+ 1 ) > 0, t i+ 1 > 0, i = 4, (3.53) 

satisfying 

Fi(ti,t i+1 ) = ti - fi(gi-i(ti),t i+1 ) = 0, t*, t*+i > 0, i = 4, (3.54) 

Therefore, in view of (13. 49p . (13. 5111 - 03.5411 . to solve the system 03.41j) 03.43^ . it is equiv¬ 
alent to solve the following single equation 


F(t n ) = t n - f n (g n -i(t, n )) = 0, t n > 0. 


(3.55) 
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In other words, to prove our proposition, we just need to show that F(-) admits a unique 
positive zero. 

We easily see that 

F(0) < 0. (3.56) 

Next we prove that 


lim F(t n ) = oo. 

t n —>oo 


(3.57) 


Noting fl3.4ip ~ fl3.43p . after a direct computation, we obtain 


lim — 

tn —^OO ~t n 

lim — 

tn—tOO tn 


tn. — 


(n ~ 1)ai2 lim K 

77/(3-11 tn—>00 t n 

1 f (i - l)(n + 2 - z)aii_i U-i 

lim 


lim 

t n —lOO 


n —1 


an i 2i(n + 1 — i) t n ->oo t. 

(i + l)(n - i)a ii+ i t i+ 1 

H- - -r— hm - 

2 i(n + 1 — i) tn-+o o t n 


hm 


dn-ln-l I 4 (n - 1) tn~>oo t n 4(n - 1) 


(3.58) 


= 2,... 

,n - 2, 

(3.59) 

na n . 

+ 

4 (n - 

-In [ 

-1)1' 

(3.60) 


Then we infer from fl3.58p ~ fl3.60p and Holder’s inequality that 


lim — < l)aii + l lim ^±4, i = 2,... ,n — 1. (3.61) 

t n —>oo t n [71 — i ~\jCLa t n —>oo t n 

In particular, we have 

lim = lim < a ™~ ln . (3. 6 2) 

tn —tOO tn —-'OO 2 _ 1 


Then using f!3.62p and Holder’s inequality again, we have 


lim 

tn~-OO 



> 

> 


1 — lim f n ^9n-l{ty)) 

tn-loo t n 

1 {F l)®r)n—1 ,. 9n—l(tn) 

1 -hm 


1 - 
1 - 


Tl&nn tn~>oo t n 

(' n ~ iKn-l 


2 T! 0 j> | n \(t nn 

n — 1 


2 n 


n + 1 
2 n 


(3.63) 


Hence (13.63p implies the desired limit f!3.57p . 

Consequently, from fl3.56p ~ fl3.57p . we conclude that F(-) admits at least one positive zero. 
In the following, we prove the uniqueness of the zero. 
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Using fl3.38j) - (l3.40j) and the constraints (I3.35P and after a direct computation, we have 


dfn—l (4—2) 4) 
dt n 


Than—In £ / J- j. \ 

4( n —1) J n—lyt'n— 2) <-nJ 


^/ P n-l(w; in-2, O 


^^n-l a n-ln-l 


(n-l) 2 A 


< 


na n —in £ / j. j. \ 

9 _i\ Jn— l\^n—2) ^n) 


4(n—1) 


9fn—l(fn—2i i 

dtn-2 


3(hi 2)ttn— In—2 J. | na n — In J- 

4(n—1) 4(n—1) 

3 (hi 2)g,yj,— In—2 


4(n 


In— In —2 £ (4 j. \ 

— 1 ) Jn —1 \^n— 2 ? ^n) 


< 


&fn—2\fn— 3? ^n—l) 

<94.-1 


/ p 2 /3T7T\ nb n -ia n -i n —1 

y 4i-n W > 4-2, 4J (n-l) 2 A 

3 (hi 2 ) 0 . 77 , — In — 2 £ /-L J- \ 

4(n—1) J hi 1 v^hi 2 5 ^7iJ 

3(71 2)0.77— In — 2 -L I nan —In -L 5 

4 ( n _l) 6 n-2 I" 4 ( n _i) 6 hi 

(hi l)o-n — 2 n— 1 


3(hl 


" 2 j" 1 /n-2(4-3,4-l) 


^4^2 (w; 4-3, 4-1) - 


4r).b„_2a rl _2n-2 


9(n—2 ) 2 A 


(yi l)ftra—2ra —1 


< 


3(n 


" 2 ^ n 1 f n -2 (4—3,4-l) 


(hi 1)0.77 —2n — 1 

3(hi—2) 

fn —2(^7i—3? ^7i—l) 


^71—1 


^71—1 


Then, from ( I3.64j) - (13.66p , we infer that 

1 (^hi) 


d fn — 1 {tn — 2 ;^n) 

0 t n 


d 4 


d fn — 1 (tn — 2 ?^n) dtn —2 
dt n —2 dt n — 1 

na n— In 


‘“n-ln 7- 7 \ 

4(n—1) 7n —1 ft-n—2jh n ; 

3(n—2)a n _p n _2 , na n —In , 

~ 4-2l ttt: TV~tn 


< 


4(n—1) 


4(n —1) 


3(n 2)a n _^77 — 2 /. /, , \ 

4(77-1) /n-l(t77-2,tn) / n _ 2 (t n _ 3 ,t n -1) 

3(n—2)a n _i n _2 + , na n —In + tn — '\ 

4 (n —1) c n—2 1 4 ( n _ 1 ) ^n 


1 - 


fn— 1(4—2)4) 4—1 9n—l(fn) 

4 4 4 

Hence, using (I3.67P and the constraints f!3.35p . we have 


dF(4) 

dl n 


= 1 

= 1 - 

> 1 - 

= 1 - 


9/„(4-i) d4- 


9t n -i dt n 

(hi 1 ) 0 . 7777—1 
hi 


£n—1 — Qn —1 (^n) 

fn (<7n —1 (4) ) d(yf„i(4) 


^/P2( w; ^_i(4)) 


4 b n an 

hiA 


dtri 


(hi 1 ) 0 , 7777—1 


fn (<7n—1 (4 ) ) Qn— 1(4) 

(n-l)o Tm - 1 ^_ i (4) tn 

fn (<?n —1 (4) ) ^(4) 


(3.64) 


(3.65) 


(3.66) 


(3.67) 


(3.68) 
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Therefore, we conclude from (J3j68j) that F(-) is strictly increasing over ( 0 ,oo), which 
implies the uniqueness of the zero of F(-). Then the proof of Proposition 13. II is complete. □ 


3.3 Constrained minimization 

By Proposition 13.11 for any w G A, we see that the equations fl3.23jl - fl3.25jl with respect to 
(ci,..., c n ) admit a solution (ci(w),..., c n (w)) given by fl3.38p - fl3.40jl . such that v defined 
by 

Vi = Wi + Ci( w), i — 1,... ,n, (3.69) 

satisfies (13.6H . 

Therefore, to find the critical points of the functional /, we consider the functional 

J(w) = l(w 1 + C\ (w), ...,w n + c n (w)) = /(w + c(w)), (3.70) 


where w 6 A. 

Noting that (13. 6 p is equivalent to (13.2211 . and multiplying (I3.22P by l r , we see that 


U r S(U- l)dx + ^ = 0. 

A 


(3.71) 


Then in view of A" 1 1 = 1, that is S' 1 P 1 1 = 1, we have 


VP- 1 (1 - U) dx + [ (U - 1 ) T S'(U - l)dx + = 0. 

Jo A 


(3.72) 


Hence from (13.7211 we may write the functional J as 


J(w) 


which is 


J(w) 


-i * r \ r i Tv* 

- Y'' / 9jW r Mdiwdx 4— / l r P _1 (1 — U) dx + b T c -, (3.73) 

2 i=1 Jn 2 J n 2 


2 „ n 

/ diW T Mdiwdx + A ^ 

i=l i=l 

n i n 


z(n + 1 — z) 
2 n 


1 _ ) dx 


(3.74) 


i=l 


i— 1 


ft is easy to see that the functional J is Frechet differentiable in the interior of A. If we 
find a minimizer w of J, which lies in the interior of A, then (w + c(w)) is a critical point 
of /. Hence, we just need to find a minimizer of J in the interior of A, denoted by int.A 
Below we first aim to find a minimizer of J in A. We begin by establishing the following 
lemma. 
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Lemma 3.1 For any w e A, there hold 


-.vV+Wi 


l,...,n, 

(3.75) 

,n. 

(3.76) 


Proof. Using (l3.38jHI3.40H and the constraints (13.3611 . we have 


e Cl < 
e Ci < 


1 (ai ( n - l)ai2 


+ 


i c 2 


an \ n 
1 


n 


a i (z - l)(n + 2 - i)au -i 

-r 3-———---e 


an \i{n + 1 — z) 2 z(rz + 1 — z) 

, (i + l)(n-z)oi i+ i c . 


•i+1 


2 z(n + 1 — i) 


, z = 2 ,..., n — 1 , 


gCn < _}_ + ( n ^-) a ^~l c C„-l 

\ 77- Ti 


Let us define an zz x n tridiagonal matrix A by 


A = 


( nan —(n — l)ai 2 
-f a 2 i 2 (n — l)a 2 2 

• • • Ajj_l 

0 
0 


0 

3(n—2) 


An 

3(n—2) 


«23 


0 

0 


\ 


-A 


zz+1 


2 In —2 2(77/ l)£7n—In—1 In 

0 (77/ l)(Z nn _l ^ 


where 


. _ (z — l)(n + 2 — z) _ ., .. 

1 = ^ I? -^zz — 7/(77/ “hi IjCLa^ 

(z + l)(n —z) 

Aj-zz+ 1 = ^ ^zz+l) 7 3, . . . , 77/ 2, 


(3.77) 


(3.78) 

(3.79) 


(3.80) 


(3.81) 

(3.82) 


and we use the notation (I3.30j) . 

We use the convention that for any two vectors a = (ai,..., a n ) T and b = (&!,..., b n y 
we write a < b if a* < 5* for all z = 1 ,..., n. 

Thus the inequalities (l3.77jMI3.79j) can be rewritten as 


A(e Cl ,..., e Cn ) r < (a 1? ..., a n ) T . 

Now we denote the adjugate matrices of A and K by 

A* = (Ah)nxn, 

IC = (AX)nxn, 


(3.83) 


(3.84) 

(3.85) 
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respectively. 

By using Holder’s inequality and an induction argument we see that all the entries of A* 
and the determinant of A are positive. Then we may express A as 


A' 1 


1 

det A 


A*, 


whose entries are all positive. Hence from (I3.83[) we obtain 


(3.86) 


( e ci \ 


/ aA 

e Ci 

< 

0j{ 

V e c -y 


) 


X 

_ j=i _ 

nan — (n—l)ai2 A£i 


X ajATj 

J = 1 


\ 


(i l)(n+2 i) A* I i(n I 1 ■ ■ A* (^+l)( n A n .. /t* 


X a j^uj 

j=i 


~(n—l)a nn -iA^_ ln +na n „A* 


(3.87) 
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Therefore, by repeatedly using Holder’s inequality, we have 


/ e Cl aA 


e 1 


\e Cn a n J 


< 


(i— 1) (n+2 — i) 


E ai%-A* 

0=1 

nanAp —(n—l)ai 2 A 21 


3 = 1 


\ 


(i+1) 2 (,t i] a ii+1 A* 


+ 1 i 


E a n ajA* ■ 
j=l 

-(n-l)a„„_iA*_ ln +na nn A* 


< 


/ n n _ \ 

1^1 FI a jj £ Klj ' 

j = l j = l 

n I 

Q ajj det K 
3 = 1 


n n _ 

1^1 fl a jj S K* 

j = 1 j = l 

™ I 

Q det K 

3 = 1 


n n _ 

|Q| PJ ajj E Knj 
j=i j=i 


V 


Q ajj det K 
3 = 1 


3 = 1 


= M 


"-i'i 


E(A- 

j=l 


v £(*-% 

v=l 


= AIA^ 1 ! = |0|1. (3.88) 


/ 


Then the proof of Lemma [3.1 1 is complete. □ 


Lemma 3.2 For any wed and (0,1), there holds 


u?+w 


*dx < 



i — 1,... ,n. 


(3.89) 


For a proof of this lemma, see mm- 

To proceed further, we need the well-known Moser-Trudinger inequality (see mm) 


e w dx < C exp 


in 



V w e fH li2 (0), 


(3.90) 


where C is a positive constant depending on fl only. 

Noting that the matrices M and S, defined by (j3.12j) and (13.151) . are both positive definite, 
we have the following coercive estimate for J. 


Lemma 3.3 For any wed there exists a positive constant C independent of X such that 

n 

A w ) > ll V ^lli - C'OnAT 1), (3.91) 

i= 1 

where ao > 0 is the smallest eigenvalue of M. 
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Proof. Since the matrices M and S, defined by (j3. 12[) and (I3.15|l . are both positive 
definite, denoting by a 0 the smallest eigenvalue of M, we have 


i =1 


z=l 


j{ w) > -y + ^bid. 

Using ^5HD - (jS35D , we have 

e Ci > 


(3.92) 


In e“i +w< dx 


2i(n + 1 - i)au 2 i(n + 1 - i) f n e 2u i +2wi dx' 

Then by the constraints (13.36|) we obtain 

2 nbj 


e Ci > 


i(n + 1 — i)\ f n e u ° +Wi dx ’ 


which implies 


2 nbi 

Ci > m —-- 

i(n + 1 — i) 


In A —In / e Ui+Wi dx , z = l,...,n. 
Jn 


(3.93) 


(3.94) 


(3.95) 


In view of Lemma 13.21 and the Moser-Trudinger inequality, we estimate the last term in 
(13.95P with 


In / e u ° +Wi dx 

Jn 

< ^ In [ e su < +sw ‘ 


1 — s 

*dx H-(In A — In4 nbi) 

n s 

s 2 \ 1- s 

< - ( In C + smaxu- + ——— 11Vxt;* 111 -f-(InA — In 4n6 ? ) 


n 


167r 


< 


1 — s 


11V it?* 11 2 H-(In A — \ndnbi) 4—In C + max iq, i — (3.96) 


Then inserting (13.961) into (13.95j) gives 

Ci > 


s „ 2 1 


Vuijllo-(hi A + InC — In 4 nbd 

167 T s 

— In 2 i(n + 1 — i) — max«-, i = 1 ,..., n. 


(3.97) 


Hence combining (I3.92j) and (13.97|) we get 


J(w) > 


do 


167T/ 


Eii v ^-tE bi (In A + In C — In 4 nbi) 


2=1 


2=1 


bi An 2i[n + 1 


i\ + max u) 


(3.98) 


2=1 


which concludes the lemma with taking s suitably small. □ 

Noting that J is weakly lower semicontinuous in A and using Lemma 13.31 we infer that 
J has a minimizer in A. 
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3.4 Interior minimizer 

In the sequel we show that the minimizer of J obtained above is an interior point of A 
when A is suitably large. To this end, we first estimate the value of the functional J on the 
boundary of A. 


Lemma 3.4 On the boundary of A there exists a constant C > 0 independent of X such that 

|fi|A 


inf J( w) > 

weM 2 


— (7(ln A + VX + 1). 


Proof. On the boundary of A, at least one of the following n conditions occurs: 

Anbun 


of = 


X 




(3.99) 


(3.100) 


Without loss of generality, if i — 1 in (I3.100p . then using (13.88j) and Holder’s inequality, 
we conclude 


e Hi < 


< 


Y a \ a jA\j 
3= 1 


nanA*^ — (n — l)ai 2 A 


21 


a 


Zj _ z u 

i n aiAPi n <*« e A 


3=2 


+ 


3=2 3=2 


2 _ 2 
an El a j 3 det K yrnri a,jj det K 

3=2 j =2 

2 n 

On v ai1 _o 

J — z 

8nb± 2(n — l)^/nbi\0,\ 


(n + 1 )A (n + 1)V\ 

Hence using Lemma [3.11 and (13.1011) . we infer that 


(3.101) 


a E 


i(n + 1 — i) 


> 


1=1 

|H|A 


2 n 


X — e *e i 


Ci„u9+Wi 


dx 


ra +1 


Anbi + [n — 1 ] y/u+JojA 


(3.102) 


For other cases, we may get similar estimates as (13.1021) . 

Then, estimating Ci,i = l,...,n, as done in Lemma 13.31 we obtain desired estimate 

dfflD. □ 

At this point, we need to find some suitable test function, which lies in the interior of 
A. We aim to compare the values of the functional at the test function with that on the 
boundary of A. 
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(3.103) 


It was proved in [5T] that for // > 0 sufficiently large, the problem 


Av = ne u °A v ( e u i +v - 1) + 


AnNi 

w 


in fl, i = 1,... ,n, 


admits a solution v?, satisfying u? + v? < 0 in fi, cf = p- f Q ifdx —> 0 , and wf = vf — cf —> 
—u® pointwise as /i —> oo, i — 1,..., n. In particular, we have the limits 


lim / e u ^ + ^dx = \Q\, lim / e 2u * +2w *dx — | 0 |, i — 

J n ^oo J n 


(3.104) 


Let us introduce an n x n tridiagonal matrix H(w M ) defined as 

/ nan — (n — 1 )|A 2 | 

-f| 0 | 2 (n — l)a 2 2 


A( w") = 


0 

,3(n-2)|^| 


0 \ 
0 


(i-l)(n+2-i) |q| _ i (i+l)(n-i) | n | 


V 


o 

o 


■^|n| 2 (n — l)a„_i n _i -f|Q| 

0 -(n-l)|Q| na nn J 


,(3.105) 


where we use the notation (I3.30P with an = an(w^),i = 1,... ,n. Then we see from (13.1041) 
that 


lim A{ w**) = \n\K. 

fl—too 


(3.106) 


It is easy to see from Jensen’s inequality that all the cofactors and the determinant of 
^(w^) are positive, which implies in particular H(w M ) is invertible and all the entries of 
J 4^ 1 (w M ) are positive. Then it follows from (I3.106P that 


lim A-\^) = —K-\ 

p-tOO 


(3.107) 


Hence we conclude from the limits (13. 1041) . (13.1071) . and the definition of A that, for a 
fixed A 0 > 0 large and for any e e (0,1), there exists a ji £ 1, such that 


= (tCi £ ,..., w^ e ) e int^4, 

(3.108) 

for every A > A 0 , and there hold 


diagjan(wf),.. .,o„ n «) j < 2|H|diag|l,..., l|, 

(3.109) 

< A-\w^) < < . 2 ,k-\ 

|H| - v ’ ~ |H| \n\ 

(3.110) 

Here we use the convention that for n x n matrices G = ( gij) n xn and H = 
G < H if gij < hij for all — 

(hij)nxn, we write 

Now we prove the following comparison result. 
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Lemma 3.5 When A > 0 is suitably large, for the test vector w Pe given by (13. 108j> . there 
holds 


J(w^) - inf J (w) < -1. 
■w&dA 


(3.111) 


Proof. For w ,,e given by (13. 108j> . using Proposition 13.11 we may get a corresponding 
vector (ci(w^),..., c n (w Me )) defined by (I3.38H — (13.401) . Then applying Jensen’s inequality 
(13.401) . we obtain 


and 


3 Cl(w'' s ) _ 


> 


Pi( w^;e C2 ( wM£ )) 
2a u (wf £ ) 

Pi(w^;e C2 ( wM£ )) 

a u (wf e ) 


1 + \ 1 — 


AbiOuiw^ 




nXPf (w^;e C2 ( wMe )) 

26, 


nAPi (w^ e ; e C2 ( w,ie )) 


Analogously, 


3 Ci(w^) > 


> 




| O | /1 n — 1 


.Me 


+ 


3 c 2 (w^) \ _ 


anjuq j \ n n 


26i 

x\n[ 


(3.112) 


_ + (i — l)(n + 2 — i) ^ Ci _,( w t*e) 


y«( w^) > 


an{wf e ) \ i(n + 1 — i) 2i(n + 1 — i) 

(z + l)(n-z) (w „ e) )_ 2 nbj 

2i(n + 1 — i) J Xi(n + 1 — i) 101 

2 b n 


M ( 1 | rt ' 1 c c- 1 ( 


a,mK £ ) v n n 


A|0| 


, i = 2,..., n-l,(3.113) 
(3.114) 


where we use the notation (13.301) with the understanding that an = au(wf e ), i — 1 ,n. 

Noting the definition (13. 105ft . we may express the inequalities (13. 1 121) — (13. 1 14ft equivalently 
as 

A( w"-) > |!1|1- 4^diag( au «),...,a„„«-)}b. (3,115) 

Since all the entries of A _1 (w Me ) are positive, we infer from (13.115ft . (13.109ft . and (13.110ft that 


(■ 


a Cl(w^ e ) C n ( W^ e 

, . . . , C 


■>)’ 


> A~‘ (w 1 *-) (|0|i - 


2 n 

x\n\ 

> A -1 (w Me ) L2| 1 —— b 


fln(wi £ ), • • • ,a nn (wlf £ ) f b 
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(3.116) 
















































Then, it follows from (13. 1 16jl that 


(l d, < |0|e+ _-^_53(/rV#. <-l.MS.1X7) 

^ ' .7—1 


At this point using (13.75H and (13.1171) we conclude that, for any small e > 0, there exists 
a positive constant C e independent of A such that 


J(w~) < (n+l)(n + 2)|n|A 

V 7 - 12 

Consequently, from (13.1181) and Lemma [3.51 we infer that 

|O|A /[n + l][n + 2] 


(3.118) 


J(w^)- inf J(w)< 

wEcM. Z 


(h±lhi±h E —l^+COnA + VX+l), (3.119) 


where C > 0 is a constant independent of A. Hence, by taking e suitably small and A 
sufficiently large in (I3.119p . we obtain the desired estimate (13.1 lip . □ 

Now we may infer from Lemma EP1 and [3~al that there exists a Ai > 0 such that, for every 
A > Ai, the functional J admits a minimizer 


w A = (tc A ,..., w^y G intA 


(3.120) 


3.5 Solution to the original system 

Since we use a constrained minimization to get a minimizer w A of J( w) = J(w + c(w)) in 
the subspace of VL 1,2 (r2), it is not obvious that whether w A gives rise to a solution of the 
system (I3.14p . Here we show that v A = (tj a , ..., u A ) defined by 

v i = w t + C;(w A ), i = l,...,n, (3.121) 

is actually a solution of the system (13.14p . 

Lemma 3.6 Let w be a minimizer of J in intA and the corresponding vector c(w) be de¬ 
termined by (I3.38p - (l3.40p . Then 


v = c(w) + w (3.122) 

must be a solution of the system (I3.14p . 

Proof. Since w is an interior minimizer of J in A, the Frechet derivative of J( w) = 
/(w + c(w)) at w should be zero, 

[d/(w + c(w))]f = 0 for any f G W h2 (Cl). (3.123) 
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By the expression of I (13. 18j) . we rewrite (13. 123j> in an explicit form 


/ fe diTMdiW + Af T US'[U - 1]\ dx 


+ [ J D f c(w)] r / ( MJS[V - 1] + dx = 0 for any f G W 1 ' 2 ^), (3.124) 


where 

D f c(w) = ^c(w + ^f)U=o, (3.125) 

is the directional derivative of c at w along the direction f, and the notation (13.5j) is used. 
Then we use (13.22j) to reduce (13.124p into 



(3.126) 


Denote by L 2 (Q) the scalar-valued or n-vector-valued function space of D-period L 2 - 
functions and decompose L 2 ( fl) as 


l 2 (D) = M n ® y, 


where 


Y= |f f GL 2 (D), J fdx = oj. 


We select a vector d G M n such that 

AILS(U - 1) + d G Y. 

Hence the relation W 1,2 (D) C Y and (I3.126P lead to 


0 = 


/ | J2 d * fT MdiW + f r (AU5[U - 1] + d) | dx 

J n |E + S) T Md t w + (f + g) r (aUS[U - 1] + d) | dx, 


for any g G M n . Consequently, we have 


J dih T MdiW + h T ^ AXIS'[U — 1] + dj j> dx = 0 for any h G Xb 1,2 (fl). 

Then we conclude from (13.1311) that w is a smooth solution of the system 

AMw = AUS(U - 1) + d, 


(3.127) 

(3.128) 

(3.129) 


(3.130) 


(3.131) 


(3.132) 
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which, after being integrated over fi, gives us 


A [\JS{U- l)dx + d|fi| = 0. (3.133) 

Jn 

Hence we infer from (13. 133|> and (13.221) that 


d = 



(3.134) 


Combining (13.132p and (13.134ft we see that v = c(w) + w is a solution of (13.14ft . Thus the 
proof of Lemma [3.61 is complete. □ 

At this stage we infer from Lemma 13.61 that when A > Ai, v A defined by (13.121ft is a 
solution of (13.14ft . Therefore part (ii) of Theorem 12.11 follows. 


3.6 Asymptotic behavior and quantized integrals 

In this subsection we study the asymptotic behavior of the solution v A of (13.14ft defined by 
(13.121ft when A —> oo and establish the quantized integrals as stated in Theorem 12.11 


Lemma 3.7 Let v A be the solution of (13.14ft given by (13.121ft . Then there holds 

lim [ (e u °i +v i — l\ dx = 0, i — (3.135) 

A->-oo V / 


Proof. Since J achieves its minimum at w A G int*4, we see from (13.118ft that, for any 
e G (0,1), there exist constants A e > 0 and C e > 0 such that 

J(w A ) = inf J (w) < ( n + 1 )(^ o +2 ) A l Q l g + C £ for all A > A e . (3.136) 

weA 12 


Noting the matrix S (defined by (13.12ft ) is positive de fini te and denoting the smallest eigen¬ 
value of S by (3q > 0, we have 


(U - 1) T S(U - l)dx > 


n p 


! U “ +Ui -1) dX- 


(3.137) 


Therefore, in view of (13.72ft . (13.74ft . and (13.137ft . and estimating c A as that in Lemma 13741 
we conclude that 


J(w A ) > 



dx — Cftln A + 1), 


(3.138) 


where C > 0 is a constant independent of A. 
Then combining (13.1361) and (13.138ft leads to 


n p 

lim sup f (- 

A->oo i=1 Jn A 


— 1) dx < 


(n + 1 )(n + 2)|0| 

6/3 0 


e for any £ > 0, 
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which implies the desired conclusion (13.1351) . The proof of Lemma [3.71 is complete. □ 
Hence part (iii) of Theorem 12.11 follows from Lemma 13.71 

To establish the quantized integrals (12.1411 . we just need to integrate the equations (13.31) 
over Q. 

The proof of Theorem 12.11 is complete. 


4 Concluding remarks 


We note that the method to establish Theorem 12.11 can be applied to prove an existence 
theorem for the problem (12.51) or (I2.6jl when the matrix K assumes a more general tridiagonal 
matrix K form, 


(l + «12 


— Oi 12 


0 \ 


K = 


1 1 + ota- 1 + cqi+1 —otu- i-i 


(4.1) 


V o 


a n n— 1 1 T C^nn—lJ 


where a 12 , a„- 1 , a ii+1) a nn - 1 > 0, i = 2 ,..., n - 1 . 

In fact, it is ready to check that all entries of K _1 are positive, A' -1 satisfies K~ l l = 1, 
and K can be decomposed as 

K = PS, (4.2) 

where P is a diagonal matrix with positive diagonal entries and A is a positive de fini te 
matrix. 

In particular, when 


«12 


Ojj—1 


Ou+ ] 


n 1 a nn —\, 

(i — l)(n + 2 — i ) 
2 


n — 1 , 


(4.3) 


the matrix K reduces to K given by (12.71) . 

Since the corresponding existence result can be stated in a similar formulation as that of 
Theorem 12.11 the details are omitted here. 

It will be of future interest to develop an existence theory when the Cartan matrix is not 
tridiagonal. 
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